Hot spots and waves in Bi 2 Sr 2 CaCu 2 08 intrinsic Josephson junction stacks - a study 
by Low Temperature Scanning Laser Microscopy 
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Recently, it has been shown that large stacks of intrinsic Josephson junctions in Bi2Sr2CaCu20g 
emit synchronous THz radiation, the synchronization presumably triggered by a cavity resonance. 
To investigate this effect we use Low Temperature Scanning Laser Microscopy to image electric 
field distributions. Apart from verifying the appearance of cavity modes at low bias we find that, 
in a high input power regime, standing-wave patterns are created through interactions with a hot 
spot, possibly pointing to a new mode of generating synchronized radiation in intrinsic Josephson 
junction stacks. 

PACS numbers: 74.50.+r, 74.72.Hs, 85.25.Cp 



Terahertz (THz) physics and technology still lacks 
good active devices (e.g. for nondestructive mate- 
rials diagnostics or chemical and biological sensing), 
a problem often referred to as the "Terahertz gap" 
Some high temperature superconductors such as 
Bi 2 Sr 2 CaCu208 (BSCCO) intrinsically form stacks of 
Josephson junctions that may operate in this regime. Re- 
cently, synchronous THz emission of almost 1000 of such 
junctions has been demonstrated 0. 

When a dc voltage V is applied to a Josephson junc- 
tion the supercurrent oscillates at a frequency / = V7$o, 
where $o is the flux quantum ("Sq 1 ~ 483.6 GHz/mV). 
Unfortunately, / is limited by the superconducting en- 
ergy gap, restricting operation of conventional Josephson 
junctions to frequencies below 700 GHz or so. Further, a 
single Josephson junction produces a small output power 
in the nanowatt or even picowatt range. Cuprate super- 
conductors, having higher energy gaps, at least in prin- 
ciple allow to operate Josephson devices up to the THz 
regime, although in reality many junction types have a 
strong internal damping, again restricting their usabil- 
ity to the sub-THz regime. The exception are intrin- 
sic Josephson junctions that are naturally formed by the 
crystal structure [1, H| . 

High frequency emission of unsynchronized intrinsic 
junctions has been observed up to 0.5 THz Q. Vari- 
ous strategies were studied to achieve synchronized THz 
radiation from intrinsic junction stacks. These studies 
included the use of shunting elements in parallel to small 
sized stacks [f| 0, H| , the excitation of Josephson plasma 
oscillations via heavy quasiparticle injection [{J hfj or the 
investigation of stimulated emission due to quantum cas- 
cade processes [III- The strategy perhaps investigated 
most intensively considered the generation of collective 
Josephson oscillations by a lattice of moving Josephson 
vortices, exciting electromagnetic cavity resonances in- 
side the stack, see e.g. [H M, M, M, M, H3, GJ, M, M, 
[2l| . Here, typically stacks or arrays of stacks with lat- 
eral dimensions of a few microns or smaller, consisting of 
some 10 junctions have been studied, perhaps with mod- 
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Figure 1: (color online). Principle of Low Temperature Scan- 
ning Laser Microscopy (LTSLM). The beam of a diode laser 
(A = 680 nm) is deflected by a scanning unit and focused onto 
the surface of the sample. In order to provide a load line for 
stable operation, the mesa (yellow structure) is biased using 
a current source and variable resistor in parallel to the mesa. 
The local heating by 0.3-3 K due to the laser causes a response 
AV serving as the contrast for the LTSLM image. In order 
to improve the signal-to-noise ratio the amplitude of the laser 
beam is modulated at a frequency 10-80 kHz and analyzed by 
a lock-in amplifier. 



est success. By contrast, in the recent experiment 
much larger structures consisting of almost 1000 junc- 
tions having lateral dimensions in the 100 /jm range have 
been used to observe coherent off-chip THz radiation 
with an output power of some /iW. Phase synchroniza- 
tion was presumably achieved via the excitation of cavity 
resonances. 

This experiment immediately triggered activities to 
understand the mechanism of THz generation in large 
intrinsic junction stacks [13, HH, HH, [25[ . Here we report 
on a study of such structures using Low Temperature 
Scanning Laser Microscopy (LTSLM), cf. FigfTJand Rcf. 
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[26| . The method allows to visualize electric field and 
current distributions with a lateral resolution of about 
1.5 /im. Standing electromagnetic waves have been im- 
aged this way. Furthermore, the large structures used 
in Q - even larger structures have been suggested [13] - 
are prone to substantial heating. LTSLM is also able to 
shine light on such heating effects. 

For the experiments Bi 2 Sr 2 CaCu208 (BSCCO) single 
crystals were grown using the floating zone technique in 
a four lamp arc-imaging furnace, and annealed at 600 °C 
and 1 atm (argon 99% and oxygen 1%). A T c of about 
83 K is typical for our samples, which are slightly under- 
doped. A gold layer 600 A thick was sputtered onto the 
surface of a cleaved BSCCO single crystal. Noticeably, 
the single crystals were cleaved in vacuum to have a 
good contact between the gold and the single crystal. 
Then conventional photolithography was used to define 
the sizes of a mesa structure in the a-b plane. The length 
of the mesas was 330 /im, the width varied between 30 /im 
and 80 /im. By using Ar ion milling, the samples were 
etched down to yield a mesa thickness of 1 /im along the 
c-axis, implying that there are about 670 intrinsic junc- 
tions involved in the stack. Polyimide was used to sur- 
round the edge of the mesa as an insulator, and a gold 
wire was attached to the mesa by silver paste. Other gold 
wires were connected to the big single crystal pedestal as 
grounds. In order to provide a load line for stable oper- 
ation, the mesas were biased using a current source and 
variable resistor in parallel to the mesa, cf. FigfTJ In 
total we measured 6 mesas on 3 different crystals. 

We first discuss data of a 40 /jm wide mesa. Figure 
[2] shows the current voltage characteristic (IVC) of this 
sample at 25 K together with LTSLM data, taken in a 
bias regime where THz emission was detected in Q • For 
voltages above 480 mV the LTSLM images are smooth, 
cf. image (G), except for an enhanced response AV along 
the edge of the mesa. Here, the laser warms up the mesa 
more effectively than in the inner region where it is at- 
tenuated by the gold layer covering the mesa surface. 
The negative sign of AV is consistent with the fact that 
the (subgap) resistance of the junctions decreases with 
increasing temperature. At lower voltages, cf. images 
(A) to (F), a wave- like structure appears which is very 
reminiscent of images of standing electromagnetic waves 
(cavity resonances) , as made by Low Temperature Scan- 
ning Electron Beam Microscop y (L TSEM) both for con- 
ventional Josephson junctions [27J and intrinsic Joseph- 
son junction stacks [2fj]. Here, in essence, the beam in- 
duced temperature rise lowers the quality factor of the 
resonance, in turn changing the voltage across the junc- 
tions. The effect is strongest when the beam heats up 
the antinodes of the standing wave and weakest when 
the beam position is at its nodes (the same physics oc- 
curs in LTSLM). Thus the period in the LTSLM image 
equals half the wave length A («120 /im in our case) of 
the cavity mode. The images were taken in a state where 
60-70% of all junctions were resistive while the others 
carried no dc voltage, as estimated from the ratio of the 



voltages on the branch investigated and the outermost 
branch. The average voltage per junction, say for im- 
age (D), is about 1 mV; using Joscphson's relation we 
estimate / w 0.5 — 0.6 THz and find a mode velocity 
c = / • A « (6 — 7) • 10 7 m/s, which is close to the highest 
possible value co/n ss 9 • 10 7 m/s, where Co is the vacuum 
speed of light and n « 3.5 is the far infrared diffraction 
index 0. Thus, the patter n imaged is consistent with 
the THz emission observed in 0, although in our case 
we see that the wave forms along the long side (in the 
fundamental cavity mode with one half-wave along the 
short side seemed to have formed). 

For this mesa we have seen similar patterns on sev- 
eral branches on the IVC (i.e. with different numbers of 
junctions in the resistive state). However, for the other 
mesas we could not find corresponding signals. This kind 
of cavity modes seems to be difficult to excite. We next 
focus on higher bias currents. Figure [3] shows data for 
a 30 /im wide mesa. At bias currents below 16 mA the 
LTSLM images, apart from edge signals, showed a broad 
response being strongest near the contact lead, cf. im- 
age (A) and the corresponding line scan. At higher bias 
currents two bright stripes appear right of the current 
lead, cf. image (B). With increasing current these two 
stripes move away from each other until, above a current 
of 22 mA, cf. image (E), the right stripe joins the edge 
signal while the left stripe continues to move towards the 
left edge of the mesa. The feature seen here is very in- 
dicative of a hot spot. For conventional superconducting 
thin films, for in-plane current flow such hot spots have 
been imaged by LTSEM [H, [H]. Here, when the beam 
touches the edge of the hot spot it increases its size, lead- 
ing to a signal AV> 0. By contrast, both in the "hot" 
and in the "cold" regions no drastic change occurs and 
AV ~ 0. In our case both in-plane and out-of-plane cur- 
rents will contribute to AV. However, we still expect the 
signal to be strongest at the edge of the hot spot, qual- 
itatively leading to the same pictures. Since for a given 
current the c-axis resistance is lower above T c than be- 
low, we expect AV < 0, as measured. Note that the IVC 
shows a small jump towards lower voltages at a bias of 
16 mA. Here, the hot spot nucleates, i.e. a small part of 
the mesa is driven to a temperature above T c . For mesas 
of larger width the hot s pot feature becomes elliptic in 
shape, cf. Fig. SI in Ref . [3l| for a 70 /im wide mesa. 

At bias currents near 23 mA, cf. images (C) to (E), 
wave-like structures appear left of the hot spot bound- 
ary. Three maxima can be seen. The distance between 
them is 55 /im. Interpreting the structures as a stand- 
ing electromagnetic wave we associate a wavelength of 
A ~ 110 /im with it. The frequency of this mode can be 
estimated from the total voltage drop across the mesa, 
V = 610 mV, yielding / w 0.44 THz and a mode velocity 
of 4.8 • 10 7 m/s, which is about a factor of 2 lower than 
cq/h. With increasing bias current the wave feature first 
disappears, but, near a bias of 28 mA, is replaced by 
another wave-like feature exhibiting 4 maxima at a dis- 
tance of 35 /im, cf. image (G). The corresponding wave 
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Figure 2: (double column): (color online) Current voltage characteristic (IVC) of a 40 x 330 jim 2 large BSCCO mesa, measured 
at 25 K, together with LTSLM data. Left graph displays the IVC on large current and voltage scales. The inset shows the 
geometry of the device. The lower graph in the middle shows an enlargement of the IVC in the region where LTSLM images 
(A) to (G) were recorded. Solid arrows denote jumps in the IVC. The upper graph on the right shows for selected bias currents, 
corresponding to images (A), (D) and (G), line scans along the long side of the mesa at half width, cf. dashed line in image 
(G). 




Figure 3: (double column): (color online) Current voltage characteristic and LTSLM data of a 30 x 330 /im 2 large BSCCO 
mesa at 50 K. Red solid arrows in the IVC denote switching processes, black arrows indicate bias points where LTSLM images 
(A) to (K) have been taken. The upper right hand graph is a 3D view of image (G) . The lower graph displays line scans along 
the long side of the mesa at half width, cf. white dashed line in image (A). Line scans are shown for images (A), (B), (F) and 
(G). 



length is A « 70 /im and the calculated mode velocity 
is 3 • 10 7 m/s. We also found standing wave patterns for 
wider mesa structures, cf. Fig. SI in Ref. [31] for a 
70 /im wide mesa. 

For an N junction stack there are N different mode 



velocities c q [H, [H, H3i the fastest of which decrease 
inversely proportional to the mode index q counting the 
number of half-waves perpendicular to the layers. Thus, 
modes with, respectively, q = 2 and 3 may have been 
excited. However, one should also consider a mode ve- 
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locity ci which is lower than cq/ti. Such a scenario is 
realistic, since even the "cold" part of the mesa is likely 
to have a temperature not too far from T c . In that 
case, the in-plane London penetration depth An can be- 
come comparable to the mesa thickness d and, as dis- 
cussed in Ref. [31], the in-phase mode velocity becomes 
Ci f=a 0.13(co/n)(d/A|i), approaching zero for T — > T c . 
For A,, w 0.5 A*m (88 % of T c for BSCCO) one obtains 
Ci « 3 • 10 7 m/s in agreement with the mode velocity of 
the resonance near 28 mA. 

Having interpreted the wave structure in terms of an 
electromagnetic wave (without having measured THz 
emission directly) we should also discuss alternative ex- 
planations. In the 70 wide mesa of Fig. SI [31] the 
wave structure appeared in a region of negative differen- 
tial resistance. One could assume that, like in a Gunn 
diode, some domains move along the junctions giving rise 
to the observed pattern. However, for the case of the 
30 fim wide mesa shown in Fig. 1 the differential resis- 
tance was positive in the regime where the wave struc- 
tures appear, making this explanation unlikely. Further, 
the resonances appeared only at bath temperatures well 
below T c , and a magnetic field of about 60 G applied par- 



allel to the layers, was sufficient to destroy the pattern, 
cf. Fig. S2 in Ref. [31]. All this strongly points to 
an interpretation in terms of the Josephson effect. Re- 
garding the role of the hot spot we note that its edge is 
typically a half wave length away from the first antin- 
ode of the wave. The hot spot seems to have an active 
role in the formation of the standing wave. As the most 
likely explanation its edge may be viewed as a resistive 
termination of the cavity formed by the cold part of the 
junction which is adjustable in space by applying vari- 
able values of bias current. It thus, in a natural way, 
combines the approaches of shunting intrinsic junctions 
and using internal cavity resonances to synchronize the 
different junctions in the stack. The effect may serve 
as an important tool to tune synchronous THz emission 
from intrinsic Josephson junction stacks. 
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Supplementary Information: 

A. Transport and LTSLM data of a 70 x 330 /im 2 

large BSCCO mesa 

Figure SI shows LTSLM and transport data for a 
70 /im wide mesa on a crystal with T c = 83 K. For cur- 
rents below 15 mA, apart from edge signals, the only 
structure is a broad spot (with AV < 0) near the cur- 
rent lead that grows in size and intensity with increasing 
current, cf. LTSLM images (A) and (B). To understand 
this signal we note that the current distributes from the 
contact lead across the mesa, with an in-plane current 
partially carried by the contacting Au layer and partially 
carried by the CuC>2 planes (as a result also the c-axis 
current density is inhomogeneous at least in the upper- 
most layers where the laser beam warms up the mesa 
most strongly). When the laser beam warms up an area 
where the supercurrent flow through the Cu02 planes is 
close to critical there will be a redistribution of currents 
and consequently a signal Ay ^ 0. By contrast, under- 
critical regions will not lead to a pronounced signal. 

For currents above 20 mA (where the current voltage 
characteristic has a kink) an elliptic spot stronger in in- 
tensity and sharper in profile forms in the center of the 
stack, cf. LTSLM image (C). This signal is the 2D analog 
of the hot spot feature discussed in the main text. For 
currents above 26 mA, a wave-like structure appears to 
the left of the hot spot, consisting of two half-waves along 
the short side of the mesa, corresponding to a wavelength 
of A y f» 70 /xm. There is also a clear modulation along the 
long side; in addition a positive signal is superimposed 
on the wave somewhat complicating the determination 
of the wavelength (we have seen such positive signals in 
several cases; the most likely explanation is that an addi- 
tional junction located close to the base crystal becomes 
resistive when the laser hits its "weakest" part). Taking 
the two bright regions in image (D) as the positions of 
the antinodes we obtain a wavelength of X x = 140 /im in 
this direction. ^From the current voltage characteristic 
we find a total voltage of 670 mV, leading to V/N w 1 mV 
and a frequency of 0.5 THz. Using / = c(A~ 2 + A~ 2 ) = 
one finds a mode velocity of about 3.1 ■ 10 7 m/s, which 
is similar to the mode velocities calculated for the more 
narrow junctions. 

B. LTSLM data of a 30 x 330 /mi 2 large BSCCO 

mesa: effect of a magnetic field 

If the observed wave pattern is associated with the 
Joscphson effect, a magnetic field applied parallel to the 
layer should affect the resonance. The field should pro- 
duce roughly one flux quantum per layer to significantly 
detune the resonance. In Figure S2 we show that a 58 G 
field, corresponding to about half a flux quantum per 
junction in the "cold" part of the mesa is able to com- 



pletely suppress the resonance. The data were taken for 
the 30 x 330 /im 2 large BSCCO mesa discussed in the 
main text. 



C. In-phase mode velocity for large junction 
numbers and temperatures close to T c 

For an N junction stack one generally finds N mode 
velocities [1, 2] c q = oj p i Xj/y/1 - 2scos{irq/(N + 1)} , 
with q = 1...N. Here, uj p i denotes the Josephson plasma 
frequency and Xj is the Josephson penetration depth; s 
is the inductive coupling parameter and can be expressed 
via 1 — 2s = (Xj/X±) 2 , where Aj_ is the out-of-plane pen- 
etration depth. When the thickness d 1 of the supercon- 
ducting layers is much smaller than the in-plane London 
penetration depth A||, for the product c = uj p iXj one 

finds c = uipiXj = (co/n)^td'/2/X\\, where t is the bar- 
rier thickness and n is the refractive index. We will fur- 
ther make use of the relation Xj = 4^ • \fsdT/2, where s 

(=1.5 nm for BSCCO) is the layer period. 

For N » 1 and (A,//A^) 2 << 1 one finds for the 
in-phase mode velocity c\. 

^ UpiXj _ UpiXj 

Cl ~ V(l - 2s) + S^JW ~ y/(\j/\±.) 2 + §ir 2 /N2 

_ u p iX,j 

v /(sd'/2A 2 | )+7r 2 /2A^ 2 

The first term under the square root dominates for 
d/X\\ > TTy/s/d', where d = Ns is the thickness of 

the mesa. In this case one finds c\ f=a (co/n)y/t]~s « 
(co/n). For BSCCO, d 1 w 0.3 nm is used often, yield- 
ing n^/s/d' w 7. Note that, since d! is not defined 
very well, this number should not be taken too liter- 
ally; if d! is closer to the interlayer distance s the ra- 
tio d/X\\ is closer to n. Using A||(T = 0) rts 170 nm, 
d must be of order 1 /im to fulfill that. If the sec- 
ond term under the square root dominates, one obtains 
ci w y/2Nc/n = {c /n)dVtd 7 /sTTX ll w 0.15(c /n)(d/A||). 
Since A|| diverges for T — > T c we see that ci approaches 
zero for T T c . For A w 0.5 fim (88% of T c for BSCCO, 
assuming a Ginzburg-Landau temperature dependence 
for the in-plane penetration depth, A|| oc 1/^/1 — T/T c ) 
one obtains ci ~ O.Scq/u, which is comparable to the 
mode velocities estimated in the main text. 
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Figure SI: Current voltage characteristic (IVC) and LTSLM data of a 70 x 330 ixm 2 large BSCCO mesa at 28 K. Solid arrows 
in the IVC denote switching processes, black dashed arrows indicate bias points where LTSLM images (A) to (G) have been 
taken. The upper right hand graph is a 3D view of image (D). The lower graph displays line scans along the long side of the 
mesa at half width, cf. white dashed line in image (A). Line scans are shown for bias points (A), (B), (C) and (F). 
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Figure S2: Suppression of the standing wave pattern by a magnetic field. The data (line scans at half width) were taken for a 
30 x 330 /im 2 large BSCCO mesa. The magnetic field was applied parallel to the short side of the mesa (see inset). 



